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Abstract 

In the present Letter we show that the concept of the generalized synchro- 
nization regime in discrete maps needs refining in the same way as it has 
been done for the flow systems [PRE, 84 (2011) 037201]. We have shown 
that, in the general case, when the relationship between state vectors of the 
interacting chaotic maps are considered, the prehistory must be taken into 
account. We extend the phase tube approach to the systems with a discrete 
time coupled both unidirectionally and mutually and analyze the essence of 
the generalized synchronization by means of this technique. Obtained results 
show that the division of the generalized synchronization into the weak and 
the strong ones also must be reconsidered. Unidirectionally coupled logistic 
maps and Henon maps coupled mutually are used as sample systems. 

Keywords: generalized synchronization, discrete maps, phase tube, 
prehistory, weak synchronization, strong synchronization 



1. Introduction 

Chaotic synchronization of nonlinear dynamical systems is an universal 
phenomenon having a large fundamental significance and different practical 
applications in all fields of science and technique jH 0, 0|. The presence of 
synchronous behavior can be observed in different mathematical, physical, 
sociological, physiological, biological and other systems. There are a lot of 
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types of chaotic synchronization such as complete, phase, generalized, noise- 
induced, lag and time scale synchronization. 

One of the most interesting types of the synchronous chaotic system 
behavior is the generalized synchronization jij. This type of chaotic syn- 
chronization is traditionally introduced for two unidirectionally coupled flow 
chaotic oscillators 0, spatially distributed media 0, Q, H| or discrete 
maps 0, [loj and means the presence of the functional relation between the 
drive and response system states. This functional relation is supposed to be 
smooth or fractal jsj , although there are no technique to find the implicit form 
of this relation (except for the complete and lag synchronization regimes). 
In the framework of the existing concept accepted generally, the strong and 
weak types of the generalized synchronization may be distinguished, accord- 
ing to the properties of the functional relation. Strong synchronization is 
assumed to correspond to the smooth map between variables of the drive 
and response systems (this regime is supposed to be observed in the case of 
complete and lag synchronization), whereas the weak one means the exis- 
tence of a fractal map between them and can be detected with the help of 
an auxiliary system approach jllj |. 

Recently, we have refined the concept of generalized synchronization in 
the flow systems and shown that the state vectors of the interacting chaotic 
systems should be considered as related with each other by the functional in- 



stead of the functional relation [12J. We have also proposed the phase tube ap- 
proach explaining the essence of generalized synchronization and allowing the 
detection of the generalized synchronization regime in many relevant physical 



circumstances including bidirectionally coupled chaotic oscillators [12 . 

Now, we have to make the next important step. The notion of generalized 
synchronization has been introduced for chaotic oscillators irrelatively of the 
type of the oscillators, it covers both the flow systems and discrete maps 



10|, [13|. At the same time, the approach proposed in [12( has been developed 
only for the flow systems. In the present Letter we extend the phase tube 
approach on discrete maps coupled both unidirectionally and mutually. As 
it would be shown bellow, in the general case the relation between states of 
the interacting discrete maps being in the generalized synchronization regime 
is analogous to the functional (as it takes place in the flow systems). As a 
consequence, the concepts of the weak and strong synchronization of chaos 
must also be reconsidered. 
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2. The theory of generalized synchronization for discrete maps 



First of all, based on the results of our previous work [12| for the flow sys- 
tems, we briefly describe the refined theory of the generalized synchronization 
regime for the discrete maps. 

The definition of the generalized synchronization regime generally ac- 
cepted hitherto is the presence of a functional relation 

y = f[x] (l) 

between the drive x and response y oscillator states 0,1^. Obviously, Eq. (JT}, 
as applied to maps, should be written in the form 

Yn = F[x„], (2) 

where x n and y n are the drive and response maps, respectively The evolution 
of the vectors x n and y n is determined by 

x„+i = H(x n , g x ), /_x 
y n+ i = G(y„, g„) + crP(x n , y„), 

where H and G are the evolution operators of the considered discrete sys- 
tems, g x and g y are the controlling parameter vectors, P denotes the coupling 
term and a is the scalar coupling parameter. Without the lack of generality 
we shall suppose below the identical dimension m of the phase space of the 
drive and response systems. 



In our work [12| we have shown for the flow systems that F[-] in Eq. (pQ) 
should be considered as a functional (contrary to a functional relation), that 
means that the system state y(t) depends not only on the state of the drive 
system x(t) but on the prehistory with the length r of the drive oscilla- 
tor x(s), t — r < s < t. From the theoretical point of view, for two coupled 
flow systems the existence of the functional relation ([T]) can be proven rig- 
orously pjj only for the unidirectional type of coupling, and, in the most 
cases, this functional relation is not continuously differentiable. As far as 
the mutual coupled flow systems are concerned, the theoretical proof men- 
tioned above becomes inapplicable. The consideration of the generalized 
synchronization of flow systems from the point of view of a functional allows 
to avoid both the uncertainty of the functional relation existence and the 
nondifferentiability feature. 
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Since the flow systems may be reduced to the discrete maps with the 
help of the of Poincare secant approach (see, e.g., [II]), the functional re- 
lation existence between system states for the generalized synchronization 
regime may be extended only for unidirectionally coupled invertible maps, 
and, again, this functional relation is fractal (i.e., it is not continuously dif- 
ferentiate) typically. As far as the the non-invertible and mutually coupled 
maps are concerned, the existence of the functional relation does not seem 
to be rigorously true. Therefore, having considered all mentioned above, 
one can come to conclusion that the prehistory should be considered in the 



same way, as it had been done for the flow systems [12|]. In terms of the 
discrete maps this circumstance may be taken into account by the following 
modification of Eq. (J2J) 

y n = F[x n ,x n _i, . . . ,x„_ x ], (4) 

where K is the discrete length of the prehistory being sufficient for the unique 
determination of the state of the response map y n . 

Let xat and be the reference points belonging to the chaotic attrac- 
tors of the drive and response maps being in the generalized synchronization 
regime, respectively. Let also 5y Jk = y 3 _ k - y N _ k and 5x Jfc = xj_ fc - x^.^ 
(k = 0, . . . , K), be the vectors characterizing the deviation of the trajecto- 
ries under consideration xj_fc, yj-k from the reference trajectories ^N-k 
and yN-k- For the neighbor point xj of the drive oscillator such that 
||£xj|| = ||5xjo|| < £ its image yj in the response system is also close to the 
reference point y^ (see {H for detail), i.e., ||<5yj|| = H^yjoll < ${ £ )- Having 
supposed that 

\\Sx Jk \\<£, k = 0,...,K (5) 
and linearized Eq. (T4]), one obtains that 

K 

Syj = ^ J XN-k F l X N, • • • , XTv^j^Xjfc, (6) 
fc=0 

where «/ Xjv _ fe is the Jacobian operator for k-th variable. Since the form of F[-] 
can not be found explicitly, Eq. OH]) may be rewritten in the form 

K 

5yj = ^2A k 5x Jk , (7) 

k=0 
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where A k = J Xjv _ fc F[x7v, • • • , Xn-k] (k = 0,...,K) are the unknown ma- 
trixes. Obviously, the coefficients of A^-matrix are determined by the whole 
set of the vectors ^n~k, ■ ■ ■ , xjv, but, since the elements of this sequence are 
uniquely connected with each other by the evolution operator (J3J), one can 
assume that depends only on xn-k, i-e., A& = Ak(x.N-K)- 

Under assumption (jSJ) made above, in view of the linearity, one can write 

<5x Jfc = B fc (x Ar _ x )5xj (8) 

[where B^^n-k) is the unknown matrix^ whose coefficients depend both 
on the reference vector ^n-k and the number k of the considered deviation 
Sx.jk], which results in 

K 

Syj = A fc (x Ar _x)B fc (x A r_x)5x J , (9) 

fc=0 

and, as a consequence, in 

5yj = C5x j7 (10) 
where C is the matrix defined as 

K 

C = ^A fe (x JV -^)B fc (x Ar -^). (11) 

fc=0 

Note, also, within the framework of the traditional concept of the gen- 
eralized synchronization implying that the states of the interacting systems 
are connected with each other by continuously differentiable functional rela- 
tion (j2J) one can obtain the similar to (1101) relationship 

5yj = C5 XJ , (12) 

with the only one difference that 

C = JF[x„,]. (13) 

Despite of the similarity of Eq. (TTUj) and Eq. ( I12p there is a great difference 
between them. Indeed, Eq. (fTQ|) has been obtained under assumption that 

Except for Bo(xjv_x) = E, where E is the identity matrix. 
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Figure 1: The dependencies of x- value of logistic map (fT5|) on the prehisory time k. The 
reference trajectory HN-k is shown by bold line, two trajectories xji_fc Xj2-k (whose end 
points x,/i xj2 are close to xjv) are shown by points • (line 1) and dashed line with points 
■ (line 2), respectively. The horizontal axis is shown in the opposite direction 



the phase trajectories x.N-k and xj_^ (k = 0, . . . , K) are close to each other 
on the whole prehisory time interval with the length K (see Eq. (JSJ), whereas 
Eq. ( IT2l) requires only the nearness of two points, x^ and xj, i.e., instead of 
Eq. (jSJ) it requires only 

||<5xj||<e. (14) 

Since for the chaotic systems the phase trajectories can both converge and 
diverge, the nearness of x^ and xj (Eq. (114jl ) does not mean the fulfilment of 
the requirement (jSJ), i.e., among the vectors xj being close to x^r only small 
part of them satisfies the requirement (j3j). This statement is illustrated in 
Fig. Q] for the logisic map 

x n+ i = ax n (l - x n ), a = 3.75. (15) 

One can see that, although both the points Xj\ and Xj2 are close to the 
reference state x^ (and for both of them requirement Eq. (JT3J) is fulfilled), 
only the point Xj\ obeys Eq. (J3J due to the nearness of the whole trajectory 
(•) to aijv-fc, whereas for the point xj 2 (■) condition (j^D fails, since its 
trajectory Xj2-k is not close to the reference one XN-k on the whole prehistory 
interval with the length K. 
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Although the coefficients of the matrixes C and C are unknown, the 
validity of both Eq. (1101) and Eq. (1121) may be verified if there are N > m 
nearest neighbors xjj of the reference vector xat and corresponding them 
vectors yjj of the response map. Note also, all vectors xjj being close to x^ 
can be used to check Eq. ffl~2l) . whereas for the examination of Eq. fflOl) only 
vectors xjj are applicable whose prehistory trajectory satisfies Eq. 

Having tested the presence of the generalized synchronization (e.g., with 
the help of the auxiliary system approach) we can pick out m nearest neigh- 
bors xjj (i = 1, . . . , m) and corresponding to them vectors yjj to determine 



the coefficients of the matrix C (or C) with the help of Eq. (fit)]) (or Eq. (fT2|) ). 
respectively) in the same way as it has been done in [l2]. Afterwards, hav- 
ing determined the coefficients of the matrix C (or C) we can now find the 
vectors Szji, (i — m + 1, . . . , N) as 

5zji = C5xjj or 5zji = C5xj; (16) 

and compare them with the vectors 5yji of the response system to check 
Eq. ffTUj) (or Eq. ( TT2"j) . To characterize the degree of closeness of the vectors 
Syji and 5zji with each other one can compute the normalized differences 

Aji " \\Sy M \\ (17) 

for each pair of vectors and build their distributions. 

So, the strategy of the investigation of the generalized synchronization 
essence may be the following. Firstly, Eq. (112j) must be checked for the set 
of vectors xjj being nearest to the reference one x^v, i.e., all points satisfying 
requirement (I14p must be used. If Eq. f fl2l) is valid, it means that in the 
generalized synchronization regime the states of discrete maps are connected 
with each other by the functional relation ([2]). Alternatively, the violation of 
Eq. (lT2"j) indicates that Eq. ([2]) being the main definition of the generalized 
synchronization concept accepted hitherto should be reconsidered. In this 
case the second step consists in the verification of Eq. (flUj) (and Eq. (j4j), 
respectively) with the help of the consideration only vectors xjj whose tra- 
jectories xj,j_fc satisfy the requirement ([5]), with the rest of the vectors xjj 
used previously to check Eq. (I12p having to be eliminated from the consid- 
eration^]. Since the length K of the prehistory (or the length of the phase 



2 For the flow system this procedure has been named as the phase tube approach 12] 
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tube) is inversely proportional to the absolute value of the largest conditional 
Lyapunov exponent X[ < 0, it may be estimated as K ~ 1/|A£|. 

In this Letter the generalized synchronization in the discrete maps is stud- 
ied for two sample systems: two unidirectionally coupled logistic maps and 
two mutually coupled Henon maps. As we will see below, the concept of 
the generalized synchronization for the discrete maps needs refining in the 
same way as it has been done for the flow systems, since, in the general case, 
for the state vectors of the interacting chaotic maps the prehistory must be 
taken into account. As a consequence, the division of generalized synchro- 
nization into weak and strong ones must also be reconsidered. At the same 
time, fortunately, this modification of the generalized synchronization con- 
cept does not discard the majority of the obtained hitherto results concerning 
generalized synchronization. 



3. Logistic maps 

As the first example we consider two unidirectionally coupled logistic 
maps: 

Vn+i = f(Vn, a y ) + or(f(x n , a x ) - f(y n , a y )), 

where f(x,a) = ax(l — x), a x = 3.75, a y = 3.79 are the control parame- 
ter values of the drive and response systems, respectively, a characterizes 
the coupling strength between systems (§]. Despite of the fact that logistic 
map is the one-dimensional discrete system, it is the etalon object of non- 
linear dynamics demonstrating the wide spectrum of interesting effects, and, 
therefore, it is typically used to study different phenomena including chaotic 
synchronization. Additionally, the logistic map belongs to the non-invertible 
discrete systems, for which the existence of the functional relation is not 
proven. Due to the one-dimensional character of interacting systems (]18|) 
the vectors in discussion given above should be replaced by scalars, whereas 
all theoretical and analytical findings remain correct. 

To detect the generalized synchronization regime we have computed con- 
ditional Lyapunov exponent for system (1181) with further refinement of the 



threshold values with the help of the auxiliary system method [111] - I n Fig- 121 
the dependence of the conditional Lyapunov exponent on the coupling pa- 
rameter u is shown. It is clearly seen that conditional Lyapunov exponent is 
negative for a G [0.12; 0.18] and a > 0.265 that is the evidence of the presence 
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of the generalized synchronization regime in these regional At that, the gen- 
eralized synchronization is close to the complete (strong) one if the coupling 
parameter is a great enough, i.e. a > 0.265, whereas for a G [0.12; 0.18] the 
detected regime corresponds to the so-called weak synchronization. It is clear 
that as in the case of the flow systems in the strong generalized synchroniza- 
tion there is no need to take prehistory into account because the drive and 
response system states are related with each other by the simple functional 
relation y n ps x n At the same time, the case of weak synchronization 
(when a G [0.12; 0.18]) demands the additional investigation. 

Without the loss of generality we fix the coupling parameter to be a = 0.14 
that corresponds to the minimum negative value of the conditional Lyapunov 
exponent (marked by arrow in Fig. |2Ja). Having assumed the value of the ac- 
curacy in Eq. (JSJ) e = 0.01 we have analyzed the influence of the length K of 
the prehistory interval on the points 5yji and normalized differences (fT7j) , 
with the reference point xn being selected randoml}0. Obviously, when 
Eq. f fTOj) (or Eq. f fT2|) ) is satisfied the distribution of normalized differences 
Aji should be the 5-function. 

Fig. El b, d,f illustrates the histograms of the normalized differences Ajj 
with the increase of the length K of the prehistory. Histograms have been 
built by A" = 1000 neighbor points being closed to each other during all pre- 
history interval with the length K. To achieve such reasonable statistics for 
a given value of accuracy e we have done L iterations the quantitative values 
of which are indicated in the caption to Fig. [21 In Fig. &\c,e,g the (x,y)- 
planes characterizing the drive and response system states for the selected 
values of the control parameters are also shown. In each figure the points 
{xji, y.ji) of the interacting systems for which requirement ([5]) is fulfilled are 
also indicated. Fig. [21 6, c corresponds to the case when the prehistory is 
not taken into account at all, i.e., K — 0. This consideration (without the 
prehistory) corresponds to the traditional concept of the generalized syn- 
chronization generally accepted hitherto. It is clearly seen that in this case 
the normalized differences Ajj are distributed uniformly over the range [0; 1] 
(Fig. 121^), at that all points in the phase space of the response system are 



3 This finding has been also verified with the help of the auxiliary system approach. 

4 It should be noted that the quantitative value of the accuracy should be a small enough 
in comparison with the amplitude of the signal from the system under study and, at the 
same time, it should be a sufficiently large to provide the reasonable statistics for a given 
time of calculation. 
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Figure 2: Dependence of the conditional Lyapunov exponent on the coupling parameter 
a (a); histograms of the normalized differences Aj; built by = 1000 neighbor points 
(b,d,f) and (x, y)-planes (c,e,g) for two unidirectionally coupled logistic maps (|18j) being 
in the generalized synchronization regime (er = 0.14, marked by arrow in this figure (a)) 
for the different lengths of the prehistory: K = (the number of iteration used to achieve 
reasonable statistics is L ~ 3.3 x 10 4 ) (b,c), K = 10 (L ~ 5.5 x 10 6 ) (d,e), K = 28 
(L ~ 6.06 x 10 10 ) (f,g). In Fig. ®c,e,g the points of the interacting systems satisfying 
requirement ([5]) are also shown 
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also allocated randomly in the wide range of the y- value variation (Fig. [21 c). 
So, we have to conclude that Eq. (lT2~j) fails and, as a consequence, the tra- 
ditional viewpoint on the generalized synchronization regime in the discrete 
systems needs refining. 

When the length of the prehistory increases, e.g., for K = 10 (Fig. [5J d, e) , 
the separate peaks in the normalized difference distribution are revealed (they 
are exist due to the inhomogeneity of the chaotic attractor), although the 
points y,ji in the phase space of the response system remain distributed in 
the wide range of the y- value as before. Finally, Fig. ^2\f,g illustrates the 
analogous distributions for the optimal length of the prehistory interval (K = 
28). In this case the distribution of the normalized differences Ajj is the 
5-function (Fig. 121/), and all considered system states (xji,yji) satisfying 
requirement (jSJ) are compressed into small neighbourhood of the reference 
point (x N ,y N ) (Fig.|2j#). 

So, having implemented the strategy developed above we can conclude 
that the relationship between states of the interacting logistic maps involves 
the prehistory of the drive system evolution in the same way as it has been 



revealed recently for the systems with continuous time (12 



4. Henon maps 

As the second example we consider two mutually coupled Henon maps. 
The mutual type of coupling between interacting systems has been selected 
for the purpose of universality, i.e. the mutual coupling is so typical as the 
unidirectional one but analysis of the generalized synchronization regime is 
performed predominantly in the unidirectionally coupled dynamical systems. 
There are also attempts to extend the concept of such phenomenon to the 



systems with a bidirectional type. For example, in our previous works [16 



17] we have shown that the generalized synchronization regime in mutually 
coupled chaotic systems could be detected by the moment of transition of 
the second positive Lyapunov exponent in the field of the negative values. 
The system under study is given by: 



x n+l ~ f{. x m x ni a x) ^{fiVni Vni a v) f( x m X ni a x))i 



X 



2 

ra+1 



bxi, 



yl +1 = f{yl, yl, a y ) + tr{f(x x n , x 2 n , a x ) - f(yl, y%, a y )), 
Vl +X = by l n , 



(19) 
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where x = (x 1 ,^ 2 ) [y = (y 1 ,?/ 2 )] are the vector-states of the first [second] 
system, f(xi,x 2 ,a) = ax\{l— Xi)+x 2 is the nonlinear function, a x = 3.16779, 



a y = 2.9, b = 0.3 are control parameters, a is the coupling parameter |!8l.ll9 



For the selected values of the control parameters generalized synchronization 
determined by the moment of the transition of the second positive Lyapunov 



exponent in the field of the negative values [16|, |20J arises at a ~ 0.035. 

Now, we fix the coupling parameter to be a = 0.2 and apply the strategy 
of the investigation of the generalized synchronization essence (see above) to 
the system under study. For the chosen value of the coupling parameter the 
weak generalized synchronization is observed in system (fl9|) . As in the case 
of the logistic maps we characterize the degree of closeness of the vectors yji 
and zji by histograms of the normalized differences (j!7p built by N = 100 
neighbor points. The quantitative values of iterations L used to achieve 
such statistics are also shown in the caption of Fig. [3j At the same time, 
contrary to the case of the logistic maps considered above, the system under 
study allows to visualize the behavior of the vectors yjj and zjj in a plane. 
Therefore, in Fig. [3] along with the distributions of the normalized differences 
A ji (Fig. [3Ja,c) the vectors yj; (o) and z^ (■) (Fig. Mb,d) of the second 
Henon map ( TT9"]) are shown. Fig. [3], a, 6 corresponds to the case when all 
the nearest vectors xjj satisfying requirement ( !T4|) with e = 0.01 are used 
(i.e., the length of the prehistory interval is K = and Eq. ( fl"2"]) is examined), 
whereas Fig.[31&,c refers to the case when the prehistory of the length K = 40 
is taken into account (in this case the requirement (E]) with e = 0.01 is 
fulfilled and Eq. fflUj) is verified). It is clearly seen that in the first case the 
normalized differences Aj, are distributed uniformly over the unit interval 
(as in the case of the logistic maps considered above) and the vectors zjj 
and yji differ from each other sufficiently that testifies the failure of the 
presence of the functional relation between the interacting system states. 
But, conversely, for the second case when the prehistory is taken into account 
the distribution of A jj is a 5-function and the calculated vectors zjj are in 
the excellent agreement with the vectors yji of the second map that confirms 
the theoretical predictions and results obtained above for the unidirectionally 
coupled logistic maps. So, for in the two-dimensional maps coupled mutually 
the vector states of the interacting chaotic systems are not also related with 
each other by the continuously different iable functional relation and, again, 
the prehistory should be be taken into account. 



12 




Figure 3: Histograms of the normalized differences Aji built by N — 100 neighbor points 
(a,c) and the vectors yji (o) and zji (■) (b,d) of the second Hennon map (|19p . a = 0.2, 
for the different lengths of the prehistory interval: K = (the number of iterations used 
to achieve reasonable statistics is L ~ 5.9 x 10 4 ) (a,6), A' = 40 (i ~ 6.2 x 10 11 ) (c,d). The 
numbers i of the vectors yjj and zj^ are shown by the italic and regular fonts, respectively 
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5. Weak and strong generalized synchronization 

In the final part of our Letter we discuss briefly the existing concept of 
the weak and strong synchronization (see, e.g., |9|) concerning the gener- 
alized synchronization regime. As it has been mentioned above, the strong 
and weak types of the generalized synchronization are typically distinguished, 
according to the properties of the functional relation between states of the 
systems. The onset of generalized synchronization is believed to be character- 
ized by an unsmooth map F that becomes smooth only at sufficiently large 
coupling strength. The synchronization types characterized by a smooth 
and an unsmooth map were called a strong and weak synchronization, re- 
spectively, with the complete and lag synchronization being a particular case 
of strong synchronization. This statement jsj was based on the calculation of 
correlation dimension (and other characteristics) of attractors in the phase 
space D © M (where © and M. are the phase spaces of the drive and response 
oscillators, respectively^. 

Indeed, if one consider the attractor of two coupled logistic maps in the 
D©R-space (see Fig. [2],c), the fractal properties of it may be easily revealed. 
At the same time, the fractality of the relationship F between states of the 
interacting systems is caused by the assumption of the existence of the simple 
function relation ([T]) between system states and neglecting the prehistory. As 
it has been discussed above, the states of the interacting systems may be not 
related with each other by the functional relation and the prehistory must be 
taken into account. To introduce the prehistory into the consideration in the 
D©R-space only vectors yjj must be used which satisfy requirement flSJ) (see 
Fig. U\g)- As one can see, in this case all considered system states (xji,yji) 
satisfying requirement ([SD are compressed into small neighbourhood of the 
reference point (xn,Un), all fractal properties disappear and the relation 
F between the drive and response system states are smooth. The same 
conclusion can be drawn not only for the logistic maps ffl8l) but for the 
general case (j3J). 

Nevertheless, the concept of the weak and strong types of the generalized 
synchronization may be used in the improved form. This improvement con- 
sists in the following. When the state of the second system y n depends on the 



5 For unsmooth map F the dimension of a strange attractor in the whole phase space 
D (BR is supposed to be larger than the dimension of driving attractor in D space, whereas 
for smooth F these two dimensions must be equal. 
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prehistory (see Eq. 01])) with the length K this type of the synchronous dy- 
namics should be considered as the weak generalized synchronization. With 
the growth of the coupling strength the required length K of the prehistory 
decreases, and, for the certain value of the coupling parameter a the length K 
becomes equal to zero and the complete synchronization regime is observed 
in the system. Since for the unidirectionally coupled oscillators the length 
of the prehistory K depends on the value of the maximum conditional Lya- 
punov exponent A^, the behavior of the prehistory length agrees well with the 
finding that the strong generalized synchronization occurs, when the maxi- 
mum conditional Lyapunov exponent drops below the minimum exponent of 



the drive system [21|. When the maximum conditional Lyapunov exponent 
becomes less then the minimum exponent of the drive oscillator, the response 
system starts to be in some sense stiff enough to follow the external signal, 
whereas the required prehistory length K is equal to zero. In this case the 
states of the interacting systems are related with each other by the contin- 
uously differentiable functional relation (J2]) that should be considered as the 
strong generalized synchronization. 

So, the division of the generalized synchronization in discrete maps on the 
weak and strong ones is certainly justified. At the same time, the difference 
between them is not determined by the type of the relation F established 
between the interacting system states (whether it is smooth or fractal), it 
is smooth in both cases, at that in the case of strong synchronization the 
interacting system states are related with each other by the functional re- 
lation (J2J), whereas in for the weak one the prehistory should be taken into 
account. 



6. Conclusions 

In conclusion, we have reported that as in the case of the flow systems 
the concept of generalized synchronization in discrete maps (coupled both 
unidirectionally and mutually) needs refining, since for the state vectors of 
the interacting chaotic systems, in general, the prehistory should be taken 
into account. We have proposed the modification of the phase tube approach 
applicable to the discrete maps and analyzed the essence of the generalized 
synchronization by means of such technique. Obtained results show that 
the division of the generalized synchronization into the weak and the strong 
ones also needs refinement, i.e. in the strong synchronization interacting 
system states are related with each other by the continuously differentiable 
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functional relation whereas in the weak one the prehistory should be taken 
into account for the analysis of the generalized synchronization regime. At 
that, both in the case of the strong and weak synchronization relation es- 
tablished between the interacting system states is smooth, and the so called 
"fractality" disappears when the appropriate consideration of the prehistory 
is made. 

At the same time, the found refinement of the generalized synchroniza- 
tion in discrete maps does not discard the majority of the obtained hitherto 
results concerning its investigation. In particular, the method of Lyapunov 
exponent computation and auxiliary system approach remains valid as before 
as well as the revealed mechanisms of the synchronous regime arising 10. loT]. 
However, this refinement has an important fundamental significance from the 
point of view of the understanding of the core mechanisms of the considered 
phenomena and is supposed to give a strong potential for new approaches 
and applications dealing with the nonlinear systems. Additionally, we expect 
that the phase tube approach gives a powerful detection and classification 
tool for the chaotic synchronization phenomenon study. 
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